clc; clear;
h = [0.5 0.25 0.125 0.1 0.05 0.025];
L2_error = [0.0241648 0.00671632 0.00189399 0.00121471 0.000304804 7.61191e-05];
plot(log(h),log(L2_error),"o-");
title("L2error order")
xlabel("log(h)");
ylabel("log(L2error)");
hold on
x = linspace(-4,-0.5);
y = 2*x;
plot(x,y);
